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We use sensitive observations of three high redshift sources; [C II] 2 P-i/2 Pi/2 fine structure 
and CO (J = 2 — > 1) rotational transitions for the z = 6.4 Quasar host galaxy (QSO) J1148+5251 
taken with the Plateau de Bure Interferometer (PdBI) and Jansky Very Large Array (JVLA) 
respectively, and [C II] and CO (J = 5 -» 4) transitions from the QSO BR1202-0725 and its 

t-H sub-millimeter companion (SMC) galaxy at z — 4.7 taken with the Atacama Large Millimetre 

Array (ALMA) and the PdBI. We use these observations to place constraints on the quantity 
Az = zco — zcr, for each source where zco and zcu are the observed redshifts of the CO 
rotational transition and [C II] fine structure transition respectively, using a combination of 
approaches; 1) By modelling the emission line profiles using 'shapelets' - a complete ortho-normal 

• ?-h set of basis functions that allow us to recreate most physical line shapes - to compare both the 

emission redshifts and the line profiles themselves, in order to make inferences about the intrinsic 
velocity differences between the molecular and atomic gas, and 2) By performing a marginalisation 
over all model parameters in order to calculate a non-parametric estimate of Az. We derive 99% 
confidence intervals for the marg inalised posterior of Az of (-1.9 ± 1.3) x 1(T 3 , (-3 ± 8) x 1CT 4 
and (-2 ±4) x 10" 3 for J1148+5251, and the BR1202-0725 QSO and SMC respectively. We show 
the [C II] and CO (J = 2 — > 1) line profiles for J1148+5251 are consistent with each other within 
the limits of the data, whilst the [C II] and CO (J = 5 — > 4) line profiles from the BR1202-0725 
QSO and SMC respectively have 65 and > 99.9% probabilities of being inconsistent, with the 
CO (J = 5 — > 4) lines ~ 30% wider than the [C II] lines. Therefore whilst the observed values of 
Az can correspond to variations in the quantity AF/F with cosmic time, where F = a 2 /fi, with a 
the fine structure constant, and fi the proton-to-electron mass ratio, of both (—3.3 ± 2.3) x 10~ 4 
for a look back time of 12.9 Gyr and of (—5 ± 15) x 10~ 5 for a look back time of 12.4 Gyr we 
propose that they are the result of the two species of gas being spatially separated as indicated by 
the inconsistencies in their line profiles. 
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I. INTRODUCTION 

One of the most important principles that underlies 
much of modern physics is that of Copernicus; that our 
location in space and time is not special, such that the 
laws of physics do not change from one point in space- 
time to another. All these laws contain within them a 
relationship to the fundamental constants; electromag- 
netic interactions involve the fine structure constant a, 
gravitational interactions depends on the gravitational 
constant G and so on. It is therefore a natural question 
to ask whether these constants are constant with cosmic 
epoch. 

This question is not a new one, being proposed some 70 
years ago by Dirac |B] and Milne [25J . Since then however 
it has become of particular interest with the development 
of physical models that seek to unify the fundamental 
forces via extra dimensions, such as Kaluza-Klein, or su- 
perstring theories (see e.g. [33] , for a review of theoretical 
models). 

These extra dimensions would have only manifested 
themselves in the first instant of time following the Big 
Bang (~ 10 -43 seconds, or energies of ~ 10 19 GeV) be- 
fore compactifying to the Planck scale during the cosmic 
expansion that followed. 

The constants of nature as we observe them in our 4 
dimensional spacetime can be shown to depend on their 
integral properties over all existing dimensions |20| . and 
as such variations in the scale factors of these now com- 
pactified extra dimensions over cosmic time could mani- 
fest themselves physically in our 4D Universe as an evo- 
lution of the fundamental constants. 

Whilst the nature of the potential evolution of these ex- 
tra dimensions is not well understand, with models exist- 
ing for increasing, decreasing, or oscillating scale lengths 
(see e.g. [9j EH [32]), a measurement showing any evolu- 
tion of the constants would provide supporting evidence 
for the existence of compact dimensions. 

The current best terrestrial estimates on the evolution 
of a and /i, the proton to electron mass ratio, come from 
Rosenband et al. [38] of \a/a\ = (1.6 ±2.3) x lCT 17 yr" 1 
and Blatt et al. [3] of = (1.6 ± 1.7) x 1CT 15 yr" 1 . 

Astrophysical methods used to acquire constraints on 
these quantities are based on comparing the differences 
between the observed line centers of absorption or emis- 
sion lines from astrophysical sources, to their expected 
values in a laboratory frame on Earth. A wide range 
of different atomic and molecular transitions have been 
used for this purpose (for a review see |10] , and for more 
recent examples from individual systems in the literature 
see Table [I]) examples of which include comparisons be- 
tween the inversion line NH3 and rotational lines from 
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CS and H2CO which resulted in the current best con- 
straints on A/i//i of —3.5 x 10~ 7 from z ~ 0.685 to 
the present [14], whilst limits on the evolution of a of 
Aa/a = (— 1.5 ± 2.6) x 10 -6 have been set using Si II 
and Fe II absorption lines from a z — 1.5 Quasar (QSO) 
PQ and of Aa/a = (— 1.7 ± 1.4) x 10~ 6 using the com- 
bination of HI and OH lines from a z — 0.765 absorp- 
tion system |13| . Evidence for a spatial variation in the 
fine structure constant has also been presented [T5J [37J , 
with the variation fitting a dipole with a significance at 
the 4.2tr level, in the direction right ascension 17.5 ± 0.9 
hours, declination —58 ± 9 degrees. 

In general, using multiple transitions from the same 
species is preferable to contending with different species, 
as long as the different transitions have differing depen- 
dencies on a, as in the former case the lines are more 
likely to arise from the same spatial location, and thus 
be less affected by random Doppler shifts of their line po- 
sitions. Such 'Doppler noise' could give rise to a velocity 
offset anywhere from ~ 10 km s _1 , comparable to the 
typical velocity dispersion of the interstellar medium in 
galaxies (e.g. [33]) up to ~ 100 km s _1 if the two species 
are located in different regions of the same galaxy. 

In [23] a method was proposed to compare the redshifts 
of far infrared (FIR) fine-structure lines and low-lying ro- 
tational transitions of the CO molecule. This combina- 
tion is sensitive to variations in the combined quantity 
F = a 2 1 'fi. 

The advantage of FIR fine-structure lines is that they 
are approximately 30 times more sensitive to variations 
in a than optical and UV transitions. In particular, 
the fine structure transition of singly ionized carbon 
2 -P3/2 P1/2 nas the combined benefit of providing 
excellent sensitivity as well as intrinsically being one of 
the brightest emission lines found in star forming galaxies 
(e.g. 0S2]). 

CO transitions then provide a second independent 
value for the redshift of the source. The frequencies of 
these emission lines are dependant on /1 such that by 
calculating the quantity Az = z rot — Zf s leads to the com- 
bined quantity F: 

Az AF 

± 1 ^avg ± 

where z avg is the weighted average of available redshift 
measurements . 

This method has been used for high redshift sources to 
place constraints on AF/F. For example in a z = 5.2 sys- 
tem Levshakov et al. [22] place limits of AF/F < 2 x 10~ 5 
using the CO (J = 7 — > 6) rotational transition combined 
with the [C I] 3 P2 — > 3 Pi fine structure transition. 

Here we apply this method to the analysis of three 
systems; Using the CO (J = 2 — > 1) rotational transi- 
tion and the [C II] 2 P3/2 ~ > 2 P\/2 fine structure transi- 
tion from the z=6.4 QSO J1148+5251 and from the QSO 
BR1 202— 0725 and its companion sub- millimetre galaxy 
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TABLE I. Existing Constraints 


on the Evolution of Physical Constants 


.Physical Quantity [A ) 


z 


A "V / V 
ZA A j A 


a 


1.84 


(5.4 ± 2.5) x 10 -6 [23] 


a 


1.15 


(-0.12 ± 1.79) x 10 -6 [30] 


a 


1.84 


(5.66 ±2.67) x 10" 6 [30] 


a 


1.58 


(-1.5 ± 2.6) x 10" 6 [1] 


/i 


Milky Way 


< 28 x 10" 9 [26] 


/' 


0.685 


(-3.5 ± 1.2) x 10 -7 [H] 


A* 


0.89 


< 9 x 10" 6 [12] 


a 2 /fi 


Milky Way 


< 3.7 x 10 -7 [25] 


a 2 /fi 


5.2 


< 2 x 10" 5 [H] 


a 2 /fi 


2.79 


(6.9 ±3.7) x 10" 6 [48] 


g p a 2 /fi 


1.17 - 1.6 


(-0.1 ±1.3) x 10 -6 [34] 


g v a 2 /ii 


1.46 (+6. 


S ± l(stat) ± 6.7 (max systematic)) xlO" 6 [16] 


r 211-85 


0.247 


(-1.18 ±0.46) x 10 -5 [H] 


r 21 1-57 

3p 1 m 


0.765 


(-5.2 ±4.3) x 10 -6 [13] 



(SMG) using the CO (J = 5 — > 4) rotational transition 
and the above [C II] line we derive values of AF/F. 

We do this firstly following a traditional approach 
where we fit a single parameterized model to both emis- 
sion lines for a given system and compare the line cen- 
troids to determine values of Az, however due to the 
very high signal to noise in some of the observations, it 
becomes increasingly possible, and important, to account 
for the non-Gaussian features in a line profile. We there- 
fore take the approach of modelling the line profiles using 
shapelets, which provide a powerful way of describing any 
physical line shape by expressing it in terms of a partic- 
ular orthonormal basis. We thus take the approach that, 
whilst the emission lines from a particular system might 
not be perfectly Gaussian, if the emission is co-located 
spatially, then they should share the same non-Gaussian 
traits. 



We also present a second, non-parametric method of 
analyzing the spectral data; performing our analysis us- 
ing the principles of Bayesian inference, we are able to not 
only quantify objectively the probability that two emis- 
sion lines from the same source share the same line shape, 
but by marginalizing over the coefficients for all possible 
models, we calculate the posterior probability distribu- 
tion for AF/F independent of any particular choice of 
model. 



In Section [IT] we describe the observations used in this 
analysis. In Section [TIT] we describe the basic mathemat- 
ical framework behind shapelets, their application to our 
data, and the marginalization process. We present our 
results from all three systems in Section [TV] and our dis- 
cussion of these results, along with our conclusions in 
Sections [V] and VI respectively. 



II. OBSERVATIONS 

The observations of the [C II] line for J1148±5251 were 
made with the Plateau de Bure Interferometer (PdB) in 
2011 [37] with a total on source integration time of 17.5 
hours. We use channels smoothed to a width of 72 MHz 
with a noise level per channel of 0.79 mJy beam -1 . An 
aperture with diameter 4 arc seconds was used to extract 
the spectra, resulting in a noise level of 1.58 mJy per 
channel and a total signal to noise ratio for the line of 
~ 15. The CO (J = 2 — > 1) line was observed with 
the Jansky Very Large Array (JVLA) in 2011. Channel 
width was set to 8 MHz, and the noise per channel of 
0.036 mJy beam -1 resulted in a total signal to noise ratio 
for the line of ~ 9. Both emission lines were imaged in a 
heliocentric velocity frame, however in order to present 
all the observations in this work in a single velocity frame 
we convert to LSR velocities by subtracting 10.27 km s -1 . 

For a detailed description of the [C II] observations 
made of the two BR1202— 0725 sources see @5]. To sum- 
marise; observations of the two lines were made with the 
Atacama Large Millimeter Array (ALMA) during 2012. 
The observations lasted 25 minutes and produced signal 
to noise ratios of ~ 80 for both lines, an order of magni- 
tude better than any previous sub-mm line or continuum 
observation of the system. The data is binned in chan- 
nels ~ 16 MHz wide for the QSO, and 32 MHz wide for 
the SMG, with noise levels per channel of 0.65 and 0.8 
mJy beam -1 per channel respectively. 

The CO (J = 5 — > 4) observations were made with the 
PdB in February 2007, for a detailed description refer to 
|39j . The channel width in both cases is ~ 30 MHz, with 
noise per channel of 0.52 mJy, resulting in signal to noise 
ratios of - 14 for both the QSO and SMG. 

The PdB observations were imaged in the LSR velocity 
frame, while the ALMA [C II] observations were made in 
a topocentric frame of reference, we therefore add 29 km 
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FIG. 1. The first 5 1-dimensional shapelet basis functions 
B n . The scaling factor is set to 200 km s _1 , with amplitudes 
normalised such that each component has a maximum at 1, 
or minimum at —1. Note the maxima/minima move further 
from the centre as n increases. 



s 1 to the [C II] observations in order to account for the 
difference between the two. 



III. SHAPELETS 

A thorough description of the Shapelet formalism can 
be found in |35j . with astronomical uses being described 
in e.g, [T7J HB |3B]- Here we give only an outline to aid 
later discussion. 

Shapelets are described by a set of dimensionless basis 
functions, which in one dimension can be written: 



(2) 



where n is a non negative integer, and H n is the Hermite 
polynomial of order n. 

These are then modified by a scale factor /3 in order to 
construct the dimensional basis functions: 



B n {x;p)=p- X ' 2 <]> n {p- X x), 



(3) 



the first 5 basis functions are shown in Fig. [T] These 
basis functions are orthonormal in that we can write 



Ax B n (x; P)B m (x; /3) = S. n 



(4) 



where 5 mn is the Kronecker delta so that we can represent 
a function f(x) as the sum: 



f( x ) = ^2a n B n (x;(3), 



(5) 



where a n are shapelet coefficients. In practice, when deal- 
ing with discretely sampled data, this feature of orthonor- 
mality depends greatly on the positions of the sample 
points, the physical extent of the sample space with re- 
spect to the scaling factor (3 and the maximum number 
of basis vectors n max used for reconstruction. For a dis- 
cussion of this see [2J. Due to the nature of our fitting 
process we can calculate the orthonormality of the set of 
basis vectors used for every model. We then define the 
quantity O max which describes the maximum deviation 
away from orthonormality for our basis vectors, i.e: 



O 



Max 



f 



N d 



aibj 



E 

(6) 

and disregard any points where O max > Othreshj where 
Othresh is a threshold value for acceptance. We find that 
for values between 1 and 10% neither the individual best 
fits nor the marginalised result vary by more than 1%. 
Below 1% a large fraction of all solutions are rejected, and 
above 10-20% the amount of power in non orthogonal 
components can result in erroneous best fits. For the 
following work we therefore set O t hresh=2%. 
The total flux F to t is given by 



o,i)£ {n} 



n 
n/2 



/qq n— even 

-°° n 

(7) 

and the r.m.s deviation of the function fix), a m is given 

by 



da; x 2 f(x) 

n— even 

nipiF- 1 23 (1 -">a„ (1 + 2n) 



n/2 



(8) 



We use this definition of the rms radius to define the 
full width at half maximum (FWHM) of a particular line 
shape as 2.3548<7 m . 



A. Model Selection 

By fitting a model to an observed spectrum using a 
set of n shapelet coefficients, we would like to be able to 
make model dependent inferences about the parameters 
of that spectrum, a task for which Bayesian Inference is 
ideal. Using the MULTINEST algorithm [Jj, we are able 
to analyse efficiently spectral data containing individual, 
or multiple emission lines, where the Bayesian evidence 
returned by the algorithm can be used to objectively per- 
form model selection between both different numbers of 
coefficients, and their values. 

We have the following scenarios: 
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1. We have a two emission lines and fit two indepen- 
dent models using the parameters i^, Pi, Ci and a 
set of rii shapelet coefficients. Here v is the central 
frequency, j3 is the scaling parameter, C is a con- 
tinuum component, n represents the set of shapelet 
coefficients used in the model and the subscript 
i = 1, 2 refers to the particular emission line. 

2. We have two emission lines and fit a joint model 
using the parameters v\, V2, /3, C, a and a set of n 
shapelet coefficients. Here v\ and V2 are the central 
frequencies of the two emission emission lines which 
are allowed to vary independently, whilst /?, the 
scaling parameter, C, a continuum component, and 
n, the set of shapelet coefficients used in the model 
are the same for both emission lines. Finally a is an 
amplitude scaling parameter representing the ratio 
of the two line peaks. 

By comparing the Bayesian evidence for these two sce- 
narios we are therefore able to objectively determine 
whether the two line shapes are consistent with one an- 
other within the limits provided by the data. In both 
these cases we also determine values of the redshift asso- 
ciated with the central frequencies of the emission lines, 
and in those cases where we fit the joint model we calcu- 
late Az — zqo — z cu, an d AF/F for the model, and the 
errors associated with these quantities. 

The only constraint that we place on the line shapes 
produced by the shapelet coefficients is that they must be 
positive for all values of velocity. We therefore reject all 
points that have any negative amplitudes by subtracting 
a large value from their log likelihood when a negative 
amplitude is present. 



B. Marginalizing over all Possible Line Shapes 

Whilst the advantages of a Bayesian approach are 
clear, in that we can robustly determine both the 
maximum number of model components supported by 
the data, and the coefficients that return the greatest 
Bayesian evidence, we are also able to marginalise over 
all the line shapes described by our shapelet coefficients 
in order to calculate the probability distribution of Az 
alone. 

The probability of a particular value of Az resulting 
from model parameters 9 is given by: 



P(Az, 9\D) oc exp [ - 0.5(D - M) T N- 1 {D - M)] , (9) 

where D is our data vector, M is the vector containing 
the model described by the parameters 9 and N is the 
covariance matrix, which in all the following work we take 
to be diagonal, with elements Na = af where Ui is the 
noise in channel i. 

The probability of a particular value of Az is therefore 
given by: 



P(Az) = J P{9\D) d9, 



(10) 



where we have integrated over all model parameters 9. In 
this way we can account for the uncertainty in choosing 
a particular model, and may also include systems where 
the line shape is not well described by a gaussian, as is 
the case for the SMG in the BR1202 system described in 
Section IVCl 



IV. RESULTS 



A. J1148+5251 



Line Profiles 



The [C II] and CO (J = 2 -> 1) spectra for J1148+5251 
are shown in Fig. [2] The evidence for the [C II] line sup- 
ports a fit with n max = 7 indicating large deviations from 
the n = term (i.e. a single Gaussian). This can be seen 
in the spectrum as a set of extended wings originating 
from strong quasar outflows 27] . The central Gaussian 
component has a FWHM of 306 ±20 km s" 1 , with an in- 
tegrated flux of 7.5 Jy km s _1 . This represents 72% of the 
total flux in the emission line, with 2.9 Jy km s _1 located 
in the non Gaussian components. The total FWHM of 
the line model is 395 ± 30 km s . We can compare this 
line width to the result of fitting only a single Gaussian 
model, which results in a FWHM of 400 ± 20 km s -1 . 
In contrast the CO (J = 2 — s- 1) line is adequately de- 
scribed by a single component model, with a FWHM of 
297 ± 49 km s _1 . This is therefore consistent within 2a 
errors with the CII line. 



2. Performing a joint fit 



As described in Section III A when performing a joint 
fit we use only the velocity range that is shared between 
both data sets and as such we only fit for the joint pro- 
file out to edge of the CO data (approximately ± 500 
km s _1 ) which means we exclude the majority of the ex- 
tended wings seen in the [C II] profile from the fit. Fig. [3] 
(left) shows the best fit joint fit profile overlaid onto both 
the [C II] and CO (J = 2 — > 1) emission lines. The con- 
sistency between the line shapes can be described quan- 
titively by comparing the Bayesian evidence for a joint 
fit where the two lines are described by the same set of 
coefficients, and the evidence when the two lines are al- 
lowed to be described by different parameters. Here the 
log evidence supports a joint profile with a difference of 
1.15, corresponding to a 76% probability that the two 
emission lines are described by the same model within 
the limits of the data. 

The joint fit is described by n max = 1 and FWHM of 
395 ± 24 km s _1 . Line centroids and redshifts are given 



() 




FIG. 2. (Left) [C II] spectrum (red) with the best fit model (magenta) overlaid and (Right) CO (J=2 — > 1) spectrum 
(blue) with the best fit model (green) overlaid for J1148+5251. In both cases zero velocity is chosen to be at the peak 
of the respective model and not at a specific redshift value. 




FIG. 3. J1148+5251: (Left) The best fit joint fit profile (magenta) overlaid onto both the [CII] (red dashed) and 
CO (J= 2 — ► 1) (blue dashed) emission lines. Amplitudes have been normalised such that the model has a peak at 1.0, 
whilst the CO (J = 2 — >■ 1) and [C II] emission lines have peaks normalised such that their height relative to the model 
is as in the best fit. Both emission lines are shown with zero velocity corresponding to their model redshift. Right: 
[C II] data (red dashed line) and CO (J = 2 — > 1) data (blue dashed line) on a redshift axis with amplitudes normalised 
as in the left panel. Redshifts have been calculated using rest frame emission frequencies of 1900.539 GHz and 230.538 
GHz for the [C II] and CO lines respectively. The solid magenta and green lines are the joint fit model line profiles at 
their best fit redshift values for the [C II] and CO (J = 2 — ¥ 1) lines respectively. 



TABLE II. J1148+5251 Emission Line Parameters from Joint 
Fit 



Line 


t'oba 


z 




GHz 




[C II] 


256.137 ± 0.010 


6.4200 ± 0.0003 


CO (J = 2 - 


+ 1) 31.081 ± 0.002 


6.4173 ± 0.0005 



in Table [TT] We find Az = (-2.7 ± 0.6 x 1CT 3 ), which, 



if the difference in redshifts were due only to changes 
in the fundamental constants and not the result of in- 
trinsic differences in the line shape would correspond to 
AF/F = (-3.6 ± 0.8) x 10" 4 , representing a 4.5cr devi- 
ation from zero. The separation in redshift is shown in 
Fig. | (Right). 

Fig [4] shows the results of marginalizing over all model 
line profiles. The left panel shows the normalised 2D 
probability distribution for n max and Az marginalized 
over all shapelet coefficients for that value of n max . The 
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J1 1 48+5251 Marginalised Posterior Distribution for A z 




-3.5 -3 -2.5 -2 

4z(10- a ) 



-□.5 



Azcicn 



FIG. 4. Left: Normalised 2D probability distribution for n max and A F/F marginalized over shapelet coefficients. 
Right: Normalised ID probability distribution for A F/F marginalized over all n max and all shapelet coefficients. 



peak is consistent with a value of Az ~ —2.4 x 10~ 3 for 
all n max with significant probabilities. Beyond this region 
the probability distribution is zero to machine precision. 
The right panel shows this 2D distribution marginalised 
over all n max . The final result is consistent with a 99% 
confidence interval for the value Az of (— 2.4±1.7) x 10~ 3 . 



B. BR1202 0725 QSO component 

1. Line Profiles 

The [C II] and CO (J = 5 ->■ 4) spectra for the BR1202 
QSO are shown in Fig. [5} The [C II] emission line is best 
modeled by a single Gaussian with a FWHM of 290 ± 4 
km s -1 , whilst the CO (J = 5 — > 4) line is described by 
a single gaussian with a FWHM of 390 ± 30 km s -1 . As 
such the two profiles are not consistent with one another 
within their errors, with the CO emission line being ~ 
30% broader than the [C II] line. 



TABLE III. BR1202-0725 QSO Emission Line Parameters 
from Joint Fit 



Line 


fobs 


2 




GHz 




[C II] 


333.723 ± 0.002 


4.69496 ± 0.00003 


CO (5 -> 


4) 101.195 ± 0.005 


4.6946 ± 0.0003 



in order to see how large an effect under fitting the data 
with a single model can have on the result. In this case 
the redshifts of the two emission lines are given by zcn— 
4.69495 ± 0.00003 and z C o= 4.6948 ± 0.0003, which 
gives us Az = (—1.4 ± 3) x 10 -4 , and correspondingly 
AF/F = (—2 ± 5) x 10~ 5 if as before we assume this 
difference in redshift is due only to an evolution in the 
fundamental constants. 

Finally in Fig. 10 we show the marginalized posterior 

-5251, from which we derive 
x 10- 4 . 



distribution Az as for J1148- 
a 99% confidence interval for Az of ( —3 ± ; 



2, Performing a joint fit 

The difference in the line widths between the [C II] and 
CO (J = 5 — > 4) emission lines is demonstrated in the dif- 
ference in the evidence values when comparing the joint 
fit to independent models, with a difference of approxi- 
mately 0.5, corresponding to a 62% probability that the 
two lines do not share the same line profile. The joint fit 
is best described by a single gaussian with FWHM 293 ± 
6 k m s~ 4 . Line centroids and redshifts are given in Table 



III 



From these we find values of Az = (—3 ± 3) x 10 4 , 
and AF/F = -5 ± 5 x 10" 5 . 

Since the evidence has a slight preference for indepen- 
dent models for the two lines, we also calculate values of 
Az and AF/F for the two independently fitted models 



C. BR1202 0725 SMG component 

1. Line Profiles 

The [C II] and CO (J = 5 -> 4) spectra for the 
BR1 202- 0725 SMG are shown in Fig. [8] Whilst the 
previous line profiles have all been roughly Gaussian in 
nature, the [C II] emission line has a much more per- 
turbed shape, with n roax of 5 required to adequately de- 
scribe it with FWHM of 700 ± 21 km s" 1 . The CO line 
however is simply described by a single Gaussian with 
FWHM 970 ±90 km s _1 , such that as with the QSO the 
CO (5 — > 4) line is approximately 30% broader than the 
[C II]. 




FIG. 5. BR1202-0725: (Left) The [C II] spectrum (red) with the best fit model (magenta) overlaid. (Right) CO (J = 
5 — > 4) spectrum (blue) with the best fit model (green) overlaid. In both cases zero velocity is chosen to be at the peak 
of the respective model and not at a specific redshift value. 
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FIG. 6. BR1202-0725: (Left) The best fit joint fit profile (magenta line) overlaid onto both the [C II] (red dotted) and 
CO (J = 5 — > 4) (blue dotted) emission lines. Amplitudes have been normalised such that the model has a peak at 
1.0, whilst the CO (J = 5 — > 4) and [C II] emission lines have peaks normalised such that their height relative to the 
model is as in the best fit. Both emission lines are shown with zero velocity corresponding to their model redshift. 
(Right) [C II] data (red dashed line) and CO (J = 5 — ¥ 4) data (blue dashed line) on a redshift axis with amplitudes 
normalised as in the left panel. Redshifts have been calculated using rest frame emission frequencies of 1900.539 GHz 
and 576.267931 GHz for the [C II] and CO lines respectively. The solid magenta and green lines are the joint fit model 
line profile at the best fit redshift values for the [C II] and CO (J = 5 — > 4) lines respectively. 



2. Performing a joint fit 

The evidence heavily favors two independent line 
shapes, with a difference in the log evidence of 20, corre- 
sponding to a probability of > 99.99% that the emission 
lines do not share the same profile. Due to the non- 
Gaussian nature of their shape, comparing the line cen- 
troids for independent models makes little sense, as the 
centers are at that stage largely arbitrary, however forc- 
ing a single profile upon the pair we find a best fit de- 
scribed by n max of 3, with FWHM of 700 ± 26 km s _1 . 
Line centroids and redshifts are given in Table |IV[ and 
correspond to values of Az — (—2.0 ± 1.2) x 10~ 3 , and 
once again assuming this change in redshift is the re- 



TABLE IV. BR1202-0725 SMG Emission Line Parameters 
from Joint Fit 

Line u ohs z 

GHz 

[C II] 333.762 ± 0.009 4.69429 ± 0.00016 
CO (5 -> 4) 101.236 ± 0.022 4.6923 ± 0.0012 



suit of only an evolution in the fundamental constants, 
AF/F = (-3.5±2.0)xKT 4 . The results of the marginal- 
ization process are shown in Fig. [lO] The posterior for 
Az displays a large number of distinct peaks resulting 
from the irregular line shapes producing a large number 
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BR1 202-0725 Marginalised Probability Distribution for A z 
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FIG. 7. Left: Normalised 2D probability distribution for n max and Az marginalized over shapelet coefficients. Right: 
Normalised fD probability distribution for Az marginalized over all n max and all shapelet coefficients. 
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FIG. 8. Left: The [C II] spectrum (red dashed) with the best fit model (solid magenta) overlaid. Right: CO (5 — > 4) 
spectrum (blue dashed) with the best fit model (solid green) overlaid. In both cases zero velocity is considered to be 
at the peak of the model. 



of likely models across a range of Az. We take the most 
Gaussian like region to represent our 99% confidence in- 
terval for Az of —1.1 ± 1.5. 



V. DISCUSSION 

A. Spatial Distributions of CO and [C II] 

When calculating the variation in F by comparing the 
line centroids of different species, one of the most impor- 
tant considerations is whether or not those two species 
have the same spatial distribution within the host galaxy. 
Recent high resolution studies (~ 50 parsec scale) of star 
forming regions in M33 [31 have shown that the CO (J 
= 2 — ¥ 1) emission line is shifted by +1.6 km s _1 relative 



to the [C II]. These authors find that significant parts of 
the [C II] are not traced by the CO because the latter 
is photo-dissociated in the low-metallicity environment 
of M33. [37] compared [C II] and CO emission on the 
scale of ~ 300 pc in the spiral arms of M31 and find that 
whilst overall they trace the same regions, their maxima 
are not coincident, with the tighter correlation being be- 
tween the [C II] and H Q , i.e. tracers of star formation. 
As such, below 300 pc we are not justified in assuming 
that the two species co-exist spatially. 

Unresolved comparisons of [C II] and CO (see e.g. [H]) 
however have shown good agreement between line profiles 
on ~ kpc scales. Our comparison of the [C II] and CO (J 
= 24 1) for J1148+5251 showed that there is a 76% 
probability that the two emission lines share the same 
profile, indicating that at least on large scales they may 
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FIG. 9. Left: The best fit joint fit profile (solid magenta) overlaid onto both the C II (red dashed) and CO (5 — > 4) (blue 
dashed) emission lines. Amplitudes have been normalised such that the model has a peak at 1.0, whilst the C05— >4 
and CII emission lines have peaks normalised such that their height relative to the model is as in the best fit. Right: 
C II data (red dotted line) and CO (5 — ¥ 4) data (blue dotted line) on a redshift axis with amplitudes normalised as in 
the left panel. The solid magenta and green lines are the joint fit model line profile at the best fit redshift values for 
the C II and CO (5 — > 4) lines respectively. 
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FIG. 10. Left: Normalised 2D probability distribution for n max and AF/F marginalized over shapelet coefficients. 
Right: Normalised ID probability distribution for AF/F marginalized over all n max and all shapelet coefficients. 



indeed trace out the same spatial volume. 

For the BR1202-0725 QSO and SMG however, we 
found -65 and > 99.9 % probabilities that the CO (J 
= 5 — > 4) and [C II] emission lines have different in- 
trinsic profiles. In both cases the CO (J = 5 —> 4) was 
approximately 30% wider than the [C II] , which suggests 
that even on large scales the spatial distribution of these 
two species is not coherent. We therefore see that with 
individual systems showing such discrepancies, any one 
estimate of Az from a single source will inherently have 
some unknown error associated with offsets in the distri- 
bution of the two species of gas. 

Several physical explanations exist to explain such 



variations; Compact starburst regions can be optically 
thick in CO, but thin in [C II], and hence the CO line 
only samples the outer regions of the galaxy leading to 
larger line widths for the CO relative to [C II]. Con- 
versely, the dust emission in such regions can become 
optically thick in the rest-frame FIR (at the frequency of 
the [C II] line), reducing it's intensity. Both [C II] and 
CO emission has been observed in the photon-dominated 
regions (PDRs) of molecular clouds [TT] . Models suggest 
a separation in the two species, with [C II] tracing the 
outer layers of molecular clouds, and CO emission com- 
ing from within. Observations of the warm interstellar 
medium in the Galaxy [IS] have also show that [C II] 
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traces a larger region than the CO. 



B. The evolution of the fundamental constants 

Table [V] summarises our results for AF/F for both 
approaches - calculating the single joint model that max- 
imises the evidence, and marginalising over all possible 
models. Errors for the single fit are given as 3er values, 
whilst the result for marginalising over the posterior is 
given as a 99% confidence interval. 

Our value of AF/F for J1148+5251 of -3.3i2.3xl0" 4 
is the result of a velocity offset with a 99% confidence 
interval of 99 ± 69 km s _1 . Thus in order to account 
for the offset from zero we would require that either 
the interstellar medium has an intrinsic velocity disper- 
sion much greater than expected, or the two species are 
positioned in such a way that their profiles are consis- 
tent, but their relative velocities are significantly offset. 
From Fig. [2] we can see that the data obtained for the 
CO (2 — > 1) line only just extends beyond the width of 
the line. There is the possibility then that the baseline 
level for this emission line has an additional component 
that is not easily visible from the data. We therefore 
allowed the continuum component of the [C II] and CO 
lines to vary independently for J1148+5251 to see how 
this affected the outcome of the fit. The result was that 
it did not; The [C II] and CO fitted continuum compo- 
nents of +0.43 ± 0.37 and +0.008 ± 0.009 respectively, 
with AF/F = (-3.3 ± 0.8) x 10~ 4 as before. 

For the BR1202-0725 QSO our value for AF/F of ( 
—5 ± 15) x 10 -5 results from a velocity offset of 15 ± 
45fcms _1 , which is consistent with our estimates of the 
intrinsic velocity dispersion of the two species. 

Finally for the BR1202-0725 SMC the best fit model 
we found a velocity offset of —120 ± 60 km s _1 , whilst 
the marginalised posterior had multiple peaks between 
a range of —375 and +30 km s -1 . This is the result of 
attempting to fit a joint model to the two emission lines 
where the evidence strongly favours different profiles for 
each. We note that the range of the offset in velocity seen 
in the marginalised posterior is roughly consistent with 
the difference in the widths of their profiles when fitted 
independently (~ 300 km s -1 ) 

If we want to reach a limit of AF/F < 10~ 6 using CO 
and [C II] we will therefore need to sample large numbers 
of high z objects in order to reduce the error associated 
with the random motions of the gas to ~ 0.3 km s" 1 
and address other factors such as differing (both gas and 
dust) opacities. To go from even a modest estimate of 
this random motion of 10 km s _1 , this would require 
~ 1000 sources, assuming Gaussian statistics. Alterna- 
tively, spatially resolved imaging could be performed on 
a number of sources, thereby providing multiple, inde- 
pendent line-of-sight measurements for each source. We 
could also choose to use two species whose spatial dis- 
tributions are believed to be more closely tied, but that 
are still detectable with high signal to noise at compa- 



rable look back times, for example WeiB et al. [IE] use 
a combination of C I and CO (7 — > 6) emission from a 
z = 2.79 QSO host galaxy to derive limits on AF/F of 
6.9 + 3.7 x 10- 6 . 

Finally if we compare our results to the proposed an- 
gular dipole distribution pH [47], J1148+251 lies in a 
region with cos6> ~ —0.9 with 9 the angular separa- 
tion with respect to the dipole, and so we would expect 

Aa/a 0.9 x 10" 5 . Given AF/F = 2Aa/a - Afi/fi, 

if we take A/j,/h = we find Aa/a = (-1.8 + 1.2) x 10~ 4 
which is consistent with the value predicted by the dipole. 
BR1 202— 0725 however lies almost orthogonal to the 
dipole with cos 8 ~ 0.2 and so does not produce any 
stringent tests on the claim. 



C. The origin of Az 

Out of the three systems investigated only the [C II] 
and CO(J = 2— >1) transition lines from J1148+5251 were 
seen to be consistent with one another within the limits 
of the data, however we note that the CO(J = 2—^1) 
emission line also had the lowest signal to noise of any 
of the data presented here. The CO (J = 5 — > 4) and 
[C II] emission lines from the BR1202-0725 QSO and 
SMG however were both seen to have intrinsically differ- 
ent profiles, where in both cases the CO (J = 5 — > 4) 
was approximately 30% wider than the [C II]. As such 
we propose that for individual systems the differences in 
the redshifts Az of the CO and [C II] emission lines can- 
not be viewed as being solely due to the evolution of the 
fundamental constants, but must include a contribution 
owing to the intrinsically different spatial distributions of 
the two species of gas within the galaxy. 



VI. CONCLUSIONS 

We have presented a series of sensitive observations of 
molecular CO and [C II] emission in 3 high redshift galax- 
ies; the QSO host J1148+5251 at z = 6.4 and the QSO 
host BR1202— 0725 at z = 4.7 along with it's companion 
SMG. We have used these to quantify the differences in 
the two emission line profiles for each system and thus put 
constraints on the quantity Az = zqo ~ z Cu m two dis- 
tinct ways; Calculating the Bayesian evidence for a joint 
model using shapelets and determining the separation in 
redshift between those models, and by marginalising over 
all model parameters to calculate the posterior distribu- 
tion associated with Az independent of any particular 
model choice. 

We found that the line profiles for the CO (5 — > 4) 
transition for the BR1202-0725 QSO and SMG were 
inconsistent with the [C II], being approximately 30% 
wider whereas the CO (2 — > 1) and [C II] emission lines 
for J1148+5251 were consistent within the limits of the 
data. This suggests that, in agreement with existing high 
resolution observations, whilst the CO and [C II] might 
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TABLE V. Final derived values for AF/F 





A Z7 1 1 Y? 

IXr 1 r 




System 


Single Joint Fit 


Marginalised Result 




± (3cr errors) 


(99% Confidence Interval) 


J1148+5251 


(-3.6 ± 2.4) x 10" 4 


(-3.3 ±2.3) x 10 -4 . 


BR1202-0725 QSO 


(-5 ±15) x 10~ 5 


( -5 ±15) x 1(T 5 


BR1202-0725 SMG 


(-4 ±6) x 10~ 4 


(-4 ±6) x 10" 3 



largely trace the same components of the host galaxy, 
there are still observable differences between their line 
profiles. As such a direct conversion from Az to the 
quantity AF/F in order to constrain the evolution of fun- 
damental constants for an individual system using these 
two species of gas is problematic, with our results being 
fundamentally limited by this unknown factor that rep- 
resents the relative distributions and velocities of the gas 
within the host galaxy. Such issues are potentially elim- 



inated either with large sample sizes to reduce the ran- 
dom error, or by observing different species that are more 
closely linked. In either case, the sensitivity of ALMA 
will allow large numbers of CO and [C II] detections at 
high redshift to mitigate the issue, allowing us to reduce 
the uncertainty in the evolution of the fundamental con- 
stants at these distant look back times to the level of 
~ 10~ 6 , on a par with methods used at lower redshifts. 
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